UNCLASSIFIED 

■n  4  2  0  6  5  8 

DEFENSE  DOCUMENTATION  CENTER 

FOR 

SCIENTIFIC  AND  TECHNICAL  INFORMATION 

CAMERON  STATION.  ALEXANDRIA.  VIRGINIA 


UNCLASSIFIED 


AIR  FORCE  SYSTEMS  COMMAND 

WRIGHT-PATTERSON  AIR  FORCE  BASE 
OHIO 


FTD-TT-63-750/1  +  2 


UNEDITED  ROUGH  DRAFT  TRANSLATION 


EFFECT  OF  ATMDSPHERIC  TURBULENCE  ON  AN  AIRPIANE 
V/ITH  FLEXIBLE  WINGS  AT  DIFFERENT  SPEEDS  OF  FLIGHT  . 

BY:  Tu.  M.  Romanovskij  and  S.  P.  Strelkov 

English  Pages :  15 

SOURCE:  Russian  Periodical,  Izvestiya  Akademll  Nauk  SSSR, 
Oideleniye  Tekhnicheskikh  Naidc,  Mekhanika  I 
Mashinostroyeniye,  Nr,  4,  1959*  pp  3  -  ID 


S/179-059-000-004 


THIS  TRANSLATION  IS  A  RENDITION  OP  THE  ORIOI* 

NAL  FOREIGN  TEXT  WITHOUT  ANY  ANALYTICAL  OR 

• 

EDITORIAL  COMMENT.  STATEMENTS  OR  THEORIES 
ADVOCATED  OR  IMPLIED  ARE  THOSE  OF  THE  SOURCE 

PREPARED  RYi 

AND  DO  NOT  NECESSARILY  REFLECT  THE  POSITION 

TRANSLATION  DIVISlOH 

OR  OPINION  OF  THE  FOREIGN  TECHNOLOGY  D|. 

FOREIGN  TECHNOLOGY  DIVISION 

VISION. 

WP.AFB,  OHia 

FTD“TT-  63-750/1  +  2 


Dat8  3  Sept  19  63 


-  EFFECT  OF  ATMOSPHERIC  TURBUIENCE  ON  AN  AIRPLANE  WITH  FIEXIBIE  WINBS 


AT  DIFFERENT  SPEEDS  OP  FUDCHT 


Tu.  M.  Romanovskiy  and  S.  P.  Strelkov 


The  varying  part  of  the  lifting  force  acting  on  an  elastic  vibrating  wing 

I 

I 

of  an  airplane  in  flight  can  be  broken  down  into  two  components :  1)  the  com¬ 


ponent  which  represents  the  varying  part  of  the  lifting  force  of  an  airplane 
wing  \mder  conditions  of  flight  in  ideally  quiet  air,  and  2)  the  random  lifting 


force  brought  about  by  atmospheric  turbulence.  If  one  limits  oneself  to  the 


theory  of  the  vibrations  of  an  elastic  airplane  in  a  turbulent  flow,  the  pzeblem 

j 

of  the  vibrations  of  an  airplane  in  a  turbulent  flow  of  air  can  be  described 


with  the  aid  of  the  following  block  diagram  (Fig.  1).  Here  the  vector  g 


characterizes  the  elastic  vibrations  of  the  airplane,  the  vector  q  deteimines 
the  aerodynamic  loads  produced  as  a  result  of  these  vibrations,  and  the  vector  p 
represents  the  aerodynamic  loads  arising  from  the  presence  of  the  turbulence 


of  the  air.  Meanwhile  all  these  values  £,  q,  and  p  are  interrelated  with  the 


aid  of  linear,  generally  integrodifferential  equations  which  depend  on  the 


horizontal  speed  of  flight. 


The  properties  of  the  airplane-air  system  change  with  the  change  in  the 


speed  of  flight.  Even  with  complete  absence  of  turbulence  at  a  determined 


value  for  the  cz*itical  speed  in  the  system,  the  block  diagram  of  idiich  is 


shown  in  Fig.  1,  there  may  appear  movements  connected  with  elastic  vibrations, 
these  movements  having  increasing  amplitude  (for  example,  wing  flutter).  Such 


instability  is  explained  by  the  fact  that  the  feedback  shown  in  the  .block 


diagram  becomes  positive  at  a  determined  value  for  the  coefficients  which 
determine  the  feedback. 


The  goal  in  this  work  is  the  investigation  of  forced  vibration  of  an 


elastic  airplane  wing  \mder  the  action  of  turbulence  as  depends  on  the  speed 
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of  flight.  The  investigation  deals  with  airplanes 
with  straight  wings  which  have  considerable  length 
in  subsonic  fli^t  speed.  Therefore  there  is  taken 
as  a  basis  of  the  computation  of  the  aerodynamic  loads 
the  theory  of  the  movement  of  a  wing  of  infinite 
length  in  a  plane  flow  [l], 

7/ith  some  limitations,  the  outside  fluctuation  action  on  the  wing  may  be 
considered  to  be  a  stationary  normal  random  process  depending  on  the  time  [2], 

By  this  same  principle,  and  because  the  whole  system  is  assumed  to  be  linear, 
one  can  make  a  description  of  its  movement  on  the  basis  of  the  correlation 
theory,  and  particularly  on  the  basis  of  the  theory  of  of  the  passing  of 
stationary  random  signals  through  linear  systems  (see,  for  example,  [3]). 

In  this  way  the  vibrations  of  the  airplane  in  the  vicinity  of  the  position 
of  equilibrium  and  the  elastic  vibrations  of  its  wings  under  the  action  of  at¬ 
mospheric  turbulence  can  be  presented  as  some  stationary,  normal,  random 
process.  The  investigation  of  this  process  amomits  to  a  determination  of 
its  different  statistical  characteristics. 

The  action  of  the  random  gusts  of  wind  on  the  wing  of  an  airplane  can 
be  presented  in  the  following  fashion.  In  its  movement  with  a  horizontal 
speed  of  U  the  airplane  goes  through  rising  and  descending  currents  of  air, 
and  the  values  for  the  velocities  of  these  currents  change  arbitrarily  both 
in  one  point  in  space  relative  to  another. and  in  time.  The  change  in  the 
value  for  the  velocity  of  a  gjst  of  wind  in  a  given  point  in  space  at  the  time 
of  an  airplane  wing’s  passing  through  it  can  be  disregarded,  Random  change 
in  the  vertical  component  of  a  turbulent  velocity  v  in  the  displacement  of  an 
airplane  from  one  point  in  space  to  another  leads  to  fluctuations  in  the  ef¬ 
fective  angle  of  attack,  and  these  determine  the  occurrence  of  random  lifting 
force.  The  influence  of  the  horizontal  component  of  the  random  velocity  of 
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gusts  of  irlnd  ordinarily  are  disregarded,  since  they  are  much  less  than  the 

speed  of  the  airplane.  Close  to  the  critical  states  the  disregarding  of  the 

i 

horizontal  component  must  be  substantiated  because  It  determines  the  parametrical 
random  action  on  the  wing  and  can  lead  to  parametrical  excitation  of  vibrations 
of  the  wing  in  the  txirbulent  flow.  As  to  the  question  of  the  action  on  the  wing 
of  the  horizontal  component  of  the  velocity,  the  authors  intend  to  devote  a 
separate  study  to  this. 

The  most  effective  action  is  brought  to  bear  on  the  airplane  by  the  large- 
scale  components  of  the  turbulence.  In  the  studies  [2,  4^  there  is  presented 
some  information  on  large-scale  turbulence.  All  the  authors  indicate  that  the 
atmospheric  turbulence  in  the  first  approximation  can  be  considered  as  isotropic. 

The  radius  of  the  correlation  of  velocities  by  various  data  fluctuates  between 

100  and  200  m.  The  figures  for  magnitudes 
of  the  spectral  densities  depend  to  a  very 
great  extent  on  the  meteorological  conditions, 
the  topography  of  the  locality  and  the  alti¬ 
tude  of  the  flight.  In  Fig.  2  there  are 
Fig.  2  presented  graphs  of  the  function  of  the 

spatial  spectral  density  s[(m2/  sec2)/(rad/m3  for  the  vertical  component  of  the 
velocity  of  the  turbulence  v,  as  depends  on  the  spatial  frequency  i^£rad/mj 
on  the  basis  of  data  presented  in  At  the  time  of  thunderstorms  the  intensity 

of  the  turbulence  may  exceed  the  indicated  values  by  one  order. 

A  nmber  of  studies  have  been  devoted  to  the  question  of  the  action  of 
atmospheric  turbulence  on  an  aircraft  using  the  methods  of  the  theory  of  random 
processes.  In  the  study  ^3']  there  is  used  the  transmissive  function  of  a 
linear  system  for  determining  the  root-mean-square  value  of  the  spectral 
characteristics  of  the  lifting  force  of  a  rigid  airplane  wing.  In  the  j 

study  it  is  shown  that  taking  into  account  the  change  in  the  velocity 

'  '  -  i 
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of  -the  turbulence  along  the  wing  leads  to  a  decrease  in  the  effect  of  the  ac¬ 
tion  on  the  wing  on  the  part  of  the  turbulence.  This  change  is  not  great, 

!  ■  ‘  ■ 

Therefore  from  here  on  it  will  be  assumed  that  along  the  wings  at  each  given 
momen  t|of  time  the  velocity  is  constant.  In  the  studies  [5*  6]  there  is  taken 

into  accovint  the  effect  of  tho  elasticity  of  the  wing  on  the  bending  in  de¬ 
termining  the  statistical  characteristics  of  normal  acceleration  of  the  bend¬ 
ing  moment  of  the  wing.  There  exist  indications  to  the  effect  that  in 

I 

the  U,  S,  A,  computations  were  made  of  the  bending  and  torsional  moments  act¬ 
ing  iJi  the  root  of  the  stabilizer  of  the  aircraft  3,s  depends  on  the  speed  of 
■flight  [7]. 

1,  Equations  of  the  System  In  writing  the  equations  which  describe  the 
vibrations  of  the  airplane  in  its  flight  the  following  scheme  was  followed. 

The  free,  airplane  flying  in  the  air  produces  elastic  vibrations  with  bending 

and  torsional  deformation  of  the  wing  with  vertical  vibrations  of  the  fuse- 

I 

lage.  In  this  scheme  the  fuselage  is  considered  to  be  a  rigid  body  with 
a  finite  mass  (comparable  with  the  mass  of  the  wing)  and  possessing  an  in¬ 
finitely  great  moment  of  inertia  relative  to  the  axis  of  rigidity  of  the  wing 
(more  precisely  a  great  moment  of  inertia  as  compared  with  the  same  magnitude 
for  the  wing).  Therefore  the  fuselage  produces  only  vertical  vibrations  as 
a  solidj  and  by  virtue  of  the  symmetry  of  the  airplane  one  may  consider  only 
one  wing  with  the  half  of  the  mass  of  the  whole  fuselage  on  one  end,  T^om 
here  on  we  will  consider  only  symmetrical  vibrations. 

As  initial  vibrations  of  the  wing,  as  is  usually  done  [8^^ there . are  taken 
the  equations  of  the  bending-torsional  vibrations  of  a  quite  long  beam 

which  practically  well  represents  the  elastic  vibrations  of  a  wing  with 

a*rrT^Z  l  ,  d*Z  d*0  , 

great  extension  fe* J  "dj*"  ^  —  A  (^>  0 +/»(*.  0 

^  f/-  f  ,  ,  d>o 

■"a;  +  +  (1.1) 

where  Z(x, t)  is  the  vertical  displacement  of  the  axis  of  rigidity  of  the 
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wling  1‘eiaCive^to  the  neutral  position,  and  0(x,t)  is  the  angle  of  turn  of  the 

I  ! 

section  of  the  wing  around  the  axis  of  rigidity.  The  positive  directions  are 

1  ' 

shown  in  Pig,  3,  All  sections  of  the  wing  are  considered  to  he  rigid  and 
nondeformable  vmder  oscillations.  The  linear  mass  of  the  wing  is  m(x)  and 

1 

the  statistical  moment  m(x)a‘(x)  and  the  moment  of  inertia  relative  to 

I  . 

the  aucis  of  rigidity  are  linear,  El  and  GI_  are  the  characteristics  of  the 
1  ^ 

rigidity  of  the  wing  to  bending  and  torsion,  respectively,  amd  f^  is  the  vary¬ 


ing  part  of  the  linear  aerodyneunio  lifting  force  and  its  moment  Mi  j  fo  and  Mo 

(  U!lt  or  rO^c.  ±  ^  ^ 

1 

are  the  random  aerodynamic  force  and  moment  which  come  into  existence  as  a 

I .  ■  .  ; 

result  of  the  turbulence, 

i 

The  system  of  equation  (l)  should  be  satisfied  vmder  the  following  con¬ 


ditions 


I 


=  0, 


(».  0 
d** 


-  ^L-0.  5».Z*+!JZ(0,<)-0  ;  (1.2) 

^where  1  is  the  length  of  the  wing  and  is  the  mass  of  the  fuselaige. 

Obtaining  a  precise  solution  of  the  system  (l,l)  appears  not  to  be  possi- 

i  '  ■ 

ble.  For  finding  approximation  solutions,  as  usual  we  will  use  the  method  of 

i 

Bubnov  and  Galerkin,  The  movement  of  the  wing  we  will  limit  by  the  following 

conditions:  the  wing  maiy  bend  symmet- 
rically  in  accordance  with  a  definite 
selected  form  and  twist  according  to 
a  definite  selected  form.  The  cvirved 
form  corresponds  to  the  form  of 
that  of  the  fundamental  of  natural 


Pig,  3 


sytneetrical  bending  oscillations  of  a  cantilever  wing,  and  the  form  of  the 
torsional  oscillations  to  that  of  the  fundamental  of  the  natural  torsional 
oscillations  of  a  wing  fixed  in  cantilever  form. 

With  these  conditions  the  solution  can  be  written 

Z(x,0  =  ao(0  +  ai(0(ZiO.  O(ic.0  =  a»(0®.(*)  (1.3) 
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^where  and  a,2  variablef*  (according  to  Galerkin).  Ueanwhile 

;z^(x)  and  ®2(*)  forms  of  the  natural  movements  of  the  wing. 

By  substituting  the  solution  in  the  form  (1.3)  in  the  equation  (l.l)  and 
using  the  known  procedure  for  the  functions  a^^  a^j  and  &2  we  will  get  the 

I  .  I  I  I 

system  ^  ^  mdx  —  ^ ^  moGjii  =  ^ /,(fx  +  ^ 

0  0  ^  0  •  ^ 

^  J  mZi*  dx  —  ^  J  moZiG,  dx  + 

0 

f  0  *  0  0“ 

'S' W"* ®** ~ ^ ~ 

i  ■  “  .  •  "  I  «  . 

-  a,  5  ^  [G/p  ^  J  %dx=  J  MA  dx+^  MA  dx. 

V  ®  .  ®  ®  . 

?o:>'  solving  the  problem  within  the  framework  of  the  correlation  theory 
it  is  sufficient  to  get  a  solution  of  the  system  (l.4)>  where  instead  of  the 
outside  random  force  f2  and  the  moment  M2  there  are  given  the  force  and  the 
moment  resulting  from  the  action  on  the  wing  of  the  airplace  of  a  disturbance 

i 

produced  by  the  presence  of  the  harmonic  wave  of  the  vertical  speed 

-  Vo  exp  [/  {(ot  -  ^)]  =  Vo  exp  [/  (  ilUt  -  ^)]  (*  =  |r)  ^  1  •  5  ) 

with  a  frequency  69  and  phase  speed  U  relative  to  the  airplane  (from  here  on 
we  will  consider  that  k  is  Strouhal’s  number). 

Since  the  vertical  speed  v  of  the  ’’cantilever  disturbance”  is  given  in  the 
form  (1.5 )j  therefore  in  the  system  studied^  as  a  result  of  its  linearity, 
there  are  possible  only  harmonic  movements.  This  means  that  in  the  case  of 
the  flight  of  an  airplane  in  ’’harmonically  disturbed"  air  there  are  formed 
stationary  harmonic  oscillations  with  a  frequency  £0*  Therefore  the  station¬ 
ary  solutions  of  the  system  (l,4)  should  be  sought  in  the  form 

flo  =  e^“'(^  +  /5o).  =  +  0,  =  ei-*  (^  + /B,)  (1.6) 

where  A^,  A2,  and  32  are  the  amplitude  coefficients  depending  on  the 

frequency  CMj  Vo,  U,  and  the  parameters  of  the  airplane. 

By  substituting  (1.6)  in  (I.4)  we  will  obtain  the  equation  for  the 
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amplitude  coefficients  in  the  fom  in  which  they  are  presented  in 

,  ■  I  ■ 

(^0  +  m  [-  ^  5  mdx  +  ^jkC[k)  «p]  +  f 

'  2  *  I  ■  .  .  ■ 

+  I'Bi)  fk  C  (k)  >tp  ^  bZidx  +  (A^  +  jB,)  [**^5  + 

•  • 

I  I  , 

+  /k  ^  rp  J  6*0,  dx(-i-  a.C  (A)]  -  npC  (A)  J  60,dr]  =  ^  K.Y  (A)  t  Mr 

•  •  • 


(1.7) 


(■^0  + /Bfi)  jkr.pC (k) -^^bZidx  +  {Ai-^  jBi)^ — A*^(mZ|*i*  + 

\  O'  f  ' 

i  ,  ®  *  , 

+  +/^i)[a*5^5  msZi%dx  —  fk^  b%Z,dx(^^  +  A,C(A))  — 

0  0 

’  r  * 

,  —  6ZiMa:]  =  itp^T(A)j6Z,dr 

®  0 
I  I 

(>lo  +  /fio)  [^7  **  5  «<»0*  +  ^  mkC  (A)  J  6*0,dr]  + 

!  '  *  »  *  I 

+  (^1  +  /^i)  A*  J  maZiO,  tix  +rp/A  A  c  (A)  A,  J  6*Z,0,dr]+ 

I  I  *  r  • 

+  (-^I  +  /^»)[— A®  5~^  /mO**<te  4-/7:  spAi  ^  6*0t*(/r  — A,C(A))  — 

”  F*  5  w]  (*)  ]  =  ’'P  ^  A,Y  (A)  (  6*0,  dx 

•  ®  0 
Here  o  is  the  density  of  the  air,  o'  is  the  averaged  distance  from  the 

center  of  the  chord  to  the  axis  of  rigidity,  C(k)  is  Theodorsen's  fxmction 

[ll#  (k)  is  Sears's  function,  and  b(x)  is  the  chord  of  the  wing  (Fig,  31)) 

^»  =?&+*)»  •-Zo'yt*  (1.7») 

By  dividing  the  imaginary  and  real  parts  in  (1.7)  it  is  possible,  in¬ 
stead  of  a  system  of  three  equations  with  complex  coefficients,  to  write  a  sys¬ 
tem  of  6  linear  algebraic  equations  rel¬ 
ative  to  the  unknowns  A^,  B^,  A-,  B^,  A^, 
and  E2.  In  the"  matrix  form  of  expression 
such  a  system  will  have  the  form 

The  coefficients  a^j^  and  b^^  are  fxmctions  of  the  munber  k.  Therefore 


Oil  •  . 

•  Oi» 

A, 

ht 

On  •  • 

-  Ot« 

B, 

bt 

0,1 . 

■  Om 

A, 

bt 

o«  .  . 

•  04I 

Bi 

* 

b* 

0*1  .  . 

•  Oj» 

At 

bt 

0*1 .  . 

•  Oh 

Bt 

bt 

(1.3) 
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also'  the  solution  of  the  sjstem  (1.8)  depends  on  k.  It  is  necessary  to  find 
Ao>  Ap  A2,  and  Bg  as  functions  of  k  or  frequency  transmission  character¬ 
istics  of  "output  values"  (a^,  aj^,  and  a^)  which  determine  the  vibrations  of 
the  aircraft  relative  to  the  "input  values"  of  the  speed  determined  by  (1.5). 

The  stationary  vibrations  of  the  airplane  in  case  of  "harmonic  distiurbance" 
to  a  great  extent  are  determined  by  the  properties  of  the  determinant  fa^^J . 

In  those  areas  of  change  of  parameters  where  this  determinant  takes  on  the 
form  of  minimum  values  there  will  generally  occur  an  increase  in  the  amplitude 
coefficients,  increase  in  the  stationary  vibrations  of  the  wing.  Such 
"resonance"  is  expressed  most  clearly  on  the  approach  to  the  conditions  of 
flutter  of  the  airplane. 

Investigation  of  the  behavior  of  the  solutions  of  the  system  (I.8)  in  the 
general  form  is  not  possible  because  of  their  extreme  bulkiness.  For  solving 
the  system  (1.8)  and  obtaining  the  values  of  the  amplitude  coefficients  for 
concrete  airplanes,  as  depends  on  the  numbers  k  and  U  and  the  parameters  of 

I 

the  airplane,  there  was  prepared  a  standard  program  by  the  aid  of  which  on  a 
rapid-action  digital  electronic  computing  machine  of  the  type  "Strela"  it 
was  possible  to  obtain  numerically  the  values  of  the  amplitude  coefficients 
with  28  dlffeient  values  for  k  and  four  values  of  U  with  =  1. 

2.  Statistical  Characteristics  of  Forced  Vibrations  of  a  Wing  lii  designing 
a  wing  for  strength  practical  interest  is  afforded  by  the  statistical  charactei^ 
istlcs  of  the  bending  and  torsional  moments  of  the  wing,  namely;  l)  the  spectral 
density  of  these  values,  2)  their  root-mean-square  values,  3)  the  number  of 
cases  of  exceeding  (by  the  indicated  values)  some  given  level  per  unit  of  time. 

For  the  spectral  densities  of  the  bending  M  and  the  torsional  *C  of  the 


moments  in  the  root  of  the  wing,  by  the  rules  of  the  correlation  theory,  we 


(2.1) 
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(2.3) 


I  vrtiere  S(k)  is 'the  spectral  densitj  of  the  vertical  canponent  of  the  velocity 

i  i  ■ 

I  of  the  tiirbulence. 

II 

.  ’  The  root-mean-square  values  of^j^and^T^  of  these  moments  are  equal  to 

'  (2.2) 

!  ■  0  .  '  • 

The  root-rmean-square  values  of  the  amplitude  of  the  bending  and  toirsional 

I  '  ■  • 

vibrations  of  the  end  of  the  wing and 27^  are  represented  in  a  similar  way 

i  S^-[5.Sz(A)dA]'‘.  (2.3] 

1  •  ■  0  •  0  . 

On  the  end  of  the  wing  (1)  =  02(1)  =  1. 

j  ■  .  ■  ‘  - 

Finally,  the  frequency  of  the  excesses  of  ^^p(^)  of  some  function  qp  (t) 

fixed  by  the  level ^  per  second  is  given  by  Rice's  formula  f3j 

i  00  00 

!  ^  e*p(- J^)[S  (k)k'-dkf‘  (2.4] 

I  ®  .  0  . 

By  a  similar  system  one  can  determine  the  statistical  characteristics  of 
i  . 

other  functions  which  describe  the  movement  of  the  airplane,  for  example,  of 

'  2  2- 
the  normal  acceleration  B^z(x,t)/^t  ,  etc. 

I 

3.  Investigation  of  the  Statistical  Characteristics  of  the  Forced  Vibrations 
of  the  Model  of  an  Airplane  For  the  investigation  of  the  statistical  character- 

1 

istlcs  with  the  aid  of  an  electronic  computing  machine  a  concrete  model  of  an 

i 

airplane  was  taken.  The  numeidcal  values  of  all  the  basic  parameters,  and  also 

I 

the  forms  of  the  bending  and  torsional  vibrations  Z3^(x)  and  02(x)  lor  the  said 

I 

model  were  taken  fromClcJ,  page  214.  The  model  is  distinguished  frtm  an  air¬ 
plane  by  the  lowered  rigidity  with  regard  to  torsion  and  bending,  and  therefore 


(2.4) 


it  has  its  own  frequency  of  bending  f^  and  torsion 


(2.4a) 


/*  =  Wt/2s  =  3.5  /#  —  S'>()/2::  =  5.043tf 

I 

The  decreased  natural  frequency  is  connected  with  the  fact  that  the  intensity 
'  of  the  t\u*bulence  spectrum  is  greatly  reduced  with  the  frequency.  It  was 

I  '* 

important  to  follow  the  effect  of  the  turbulence  in  the  case  where  the  frequency 
of  the  flutter  of  the  model  still  lies  in  that  area  of  frequencies  of  the 
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Fig.  4 

Half  of  mass  of  airplane 
model  M  =  2063  kg 
Statistical  moment  of  the 
wing  S*  =  77  kgm 
Moment  of  inertial  of  the 
ydjig  3^0  =  kgm^ 


turbulence  spectra  where  the  intensity  is 
relatively  great. 

Let  us  present  some  numerical  values 
of  the  basic  coefficients  of  equation  (1.7) 
of  the  model  under  investigation. 

Average  distance  from  rigidity  line 
to  central  line  a  =  2tr/h^  ~  -0,46 
Average  statistical  moment  in  the 
end  of  the  wing  Mo  “  6,000  kOu 
Half  of  wing  span  1  =  6.12  m 


Average  chord  =  2.3  ni 

The  magnitudes  of  the  xmited  masses  of  air  and  their  inertia  are  included 
in  the  respective  coefficients. 

There  were  computed  the  squares  of  the  modules  of  the  transmission  character¬ 
istics  for  the  bending  and  torsion  of  the  wing 

+/■■»>.  (2.4b) 


of  the  model  under  investigation  with  several  values  for  the  speed  (the  computa¬ 
tion  was  done  on  the  "Strela”  in  the  Computation  Center  at  the  Moscow  State 
University;  for  obtaining  108  solutions  of  the  system  of  equations  of  the  6th 
order  (13)  five  minutes  were  required). 
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In  order  to  determine  the  influence 
of  torsion  on  these  transmission  char¬ 
acteristics  analogous  computations  were 
made  for  the  same  wing  in  the  proposition 
e2(x)  HO.  In  Fig.  4  there  are 
presented  graphs  of  the  function  ®]^(k) 
for  the  system  with  6  degrees  of 
freedom.  In  Fig.  5  there  are  the 
corresponding  graphs  of  the  system  in 
the  piroposition  02  H  0  (four  degrees 
of  freedom).  The  critical  rate  of  the 
flutter  U^,  computed  for  the  model, 
amounted  to  about  m/sec. 

From  a  comparison  of  the  trans¬ 
mission  functions  (Figs.  4  and  5)  and 
the  curves  of  the  spectral  densities 
of  the  bending  of  the  end  of  the  wing 
for  the  systems  with  six  and  four 


degrees  of  freedom  (Figs.  6  and  7)  constructed  for  different  values  of  U  m/sec, 

I 

I 

it  follows  that  taking  into  account  the  torsion  of  the  wing  leads  to  a  sharp 
accentuation  of  the  curves  of  the  transmission  functions,  and,  accordingly, 
of  the  curves  of  the  spectral  densities  with  an  increase  in  the  rate  of  U. 

In  the  tables  1  and  2  there  are  presented  confrontations  of  statistical 
characteristics  obtained  on  the  basis  of  the  curves  of  spectral  densities  by 
the  formulas  (2.2),  (2.3)*  and  (2.4)  for  the  case  =  117  m/sec.  With  the 
aid  of  the  data  of  these  tables  one  can  follow  the  change  in  the  jroot-mean- 
square  of  the  values  andSg,  and  of  the  f^mctions  N2(^)  and  Ng(^) 
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with  the  change  in  the  rate  of  U.  In 


I 

Displacement  —  bending  (^is  the  deflec—  both  the  cases  considez^d  with  an 


tion  of  the  end  of  the  wing  in  cm) 


u 

M 

ceic 

CM  • 

. 

kTm 

J»,«) 

e-0 

t-»i 

50 

1.18 

342 

1.4 

0.45X10-* 

0.28xlO-*» 

80 

1.61 

487 

1.8 

0.81  Xl0-‘ 

0.54X10-“ 

115 

2.54 

735 

1.6 

0.47 

0.47xl0-» 

(1) 

0 

1160 

>4350  ^ 

increase  in  U  there  occurs  an  Increase 
in  the  computed  root^ean-square  values 
and  the  numbers  of  intersections  of  a 


Key:  (l)  moment  in  the  root^^  the 
wing  with  a  given  bending^  of  the 
end  of  the  wing  (kg-m) 


the  root-mean-square  values  of 
with  the  increase  in  the  rate  of  U 
from  50  to  115  m/sec  are  Increased 
by  a  factor  of  12,  while  in  the  second 
case  (displacement  -  bending)  is 
increased  only  by  a  factor  of  2  with 


the  same  change  in  the  speed.  The  nimiber  of  cases  of  exceeding  a  given  level 
with  the  growth  in  speed  changes  still  more, 

T  a  b  1  e  2 

Displacement  -  bending  -  torsion  is  the  deflection  of  the  end  of  the 
wing  in  cm,  ^  is  the  angle  of  turn  of  the  wing  at  the  end) 
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4.5 

4.35 

2.02 
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44 

220 

Key;  (l)  degree;  (2)  degrees;  (3)  moment  in  the  root  of  a  wing  section 
with  a  given  angle  of  torsion^^^  on  the  end  (kgm) 

Analogous  computations  were  made  for  the  same  model  of  airplane  with  changed 

_  .  ,  ,  7 

natural  frequency  of  torsion  (f^  =  3*5  cps,  f^,  =  6.0  cps,'  U*  =  200  m/sec). 
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Meanwhile  the  character  of  the  change  of  the  statistical  characteristics  proved 

i 

to  be  about  the  same  as  in  the  first  case. 

\ 

Conclusions  The  described  method  makes  it  possible  to  obtain  simpl7  the 
statistical  characteristics  of  forced  vibrations  of  an  elastic  airplane  wing 
under  the  action  of  atmospheric  turbulence.  With  its  help  one  can  take  into 
account  vibrations  with  different  degrees  of  freedom.  The  increase  in  the 
number  of  degrees  of  freedom  per  unit  leads  to  an  increase  in  the  order  of 
solvable  system  of  algebraic  equations  by  two  units.  Meanvdiile  into  the 
standard  program  it  is  necessary  to  introduce  insignificant  changes;  the  time 
for  the  computation  correspondingly  is  increased  by  a  small  factor. 

Comparison  of  the  statistical  characteristics  of  vibrations  of  one  and  the 
same  model  of  airplane  shows  that  with  sufficiently  low  frequencies  of  torsion 
of  a  wing  close  to  the  critical  rates  of  flutter  one  should  not  disregard  the 
torsion  of  the  wing,  but  it  is  necessary  to  take  it  into  consideration  in 
computing  the  Joint  bending-torsional  vibrations  of  the  wing. 

The  authors  express  their  deep  gratitude  to  V.  A.  Dnizhinnaya  and  V.  B, 
Glasko  for  setting  up  a  standard  program  and  carrying  out  the  computations  on 
the  "Strela”  machine. 
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